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Recent progress in lattice calculations of properties of open-charm 
mesons, both regular and exotic, is reviewed, with an emphasis on spec¬ 
troscopy. After reviewing recent calculations of excited state energy lev¬ 
els I will discuss progress in extracting hadronic masses and widths of 
charmed states from Lattice QCD simulations including low-lying scat¬ 
tering channels directly, to determine phase shift data and bound state/ 
resonance properties. With regard to other properties results from recent 
calculations of the DD*ti and DDp, D*D*p couplings are presented. Be¬ 
yond regular mesons, searches for explicitly exotic (tetraquark) states are 
also reviewed. 
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1 Introduction 


Lattice QCD, the simulation of Quantum Chromodynamics (QCD) on a space-time 
grid, is commonly regarded as the ideal tool to investigate the non-perturbative 
physics of the strong interaction at hadronic energy scales. Indeed there has been 
impressive progress in recent years, with simulations at or close to physical quark 
masses with dynamical light, strange and even charm quarks. For ground state ob¬ 
servables involving mesons only, many quantities have been calculated with both full 
control of the relevant systematic uncertainties and with an impressive statistical and 
systematic precision. The masses of a number of low-lying hadrons have been de¬ 
termined in this way BEI including ground states with charmed quarks. Example 
of quantities determined with full control of systematic uncertainties are collected in 
the FLAG review [3]. The review collects lattice results relevant for flavor physics 
and tries to provide the current best lattice values for particular quantities listing the 
various lattice results categorized by the number of dynamic flavors and rating them 
by fixed quality criteria. 

For properties of hadronic excitations or for hadrons close to multi-hadron thresh¬ 
olds the situation is however complicated, and fully systematic determinations of 
excited state properties are still lacking. The reason for this is that excited states are 
much more challenging, both conceptually and computationally: 

• Looking at Euclidean-space correlation functions, excited state contributions 
appear as subleading exponentials 

(|O 2 (t)Oi( 0 )^ oc < 0 |O 2 |n >< n|Oi|0 > , 

n 

where En is the energy of the n-th state. For a given channel the sum in principle 
consists of all hadronic (including multi-hadron) states with the same quantum 
numbers. This makes the reliable extraction of excited state energy levels in 
Monte-Carlo simulations with non-vanishing statistical uncertainty very chal¬ 
lenging. 

• For hadronic resonances and close-to-threshold bound states, scattering phase- 
shifts and consequently resonance and bound state poles are not directly acces¬ 
sible in Euclidean space |1], but instead have to be extracted using the finite 
volume method pioneered by Liischer mum- 

• In determining energy levels from lattice calculations using quark-antiquark 
interpolating helds for mesons and 3-quark interpolating helds for baryons, 
it became obvious [HiEiEniin] that a basis of only quark-antiquark or only 
scattering operators is not very useful and is probably not suitable to obtain 
the true energy levels. Modern calculations therefore employ a mixed basis 
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Figure 1: D and Dg meson spectra by the Hadron Spectrum Collaboration [12]. The 
quantum numbers are denoted on the horizontal axis. All results are for up and 
down quark masses corresponding to 400MeV pions. Lattice data in green compared 
to mesons known from experiment in black. 


of both scattering and single-hadron interpolating helds, which for many in¬ 
teresting quantum numbers necessitates the inclusion of quark-line diagrams 
with backtracking quark lines, (inspired by single hadrons) often referred to as 
“disconnected diagrams”. These are especially challenging numerically. 

Consequently most calculations of charmed excitations are of an exploratory na¬ 
ture and currently lack a detailed estimate of the systematic uncertainty in the calcu¬ 
lation. A number of results reported in this review are also preliminary in the sense 
that they have not yet been published. 


2 Single hadron spectra 

As outlined in the previous section, one of the challenges for lattice calculations of ex¬ 
cited state masses consists of reliably determining excited state energy levels. Figure 
[T] shows the energy levels determined by the Hadron Spectrum Collaboration from 
single meson interpolating helds for charmed and charmed-strange mesons [T2|. The 
impressive results demonstrate that a large number of energy levels in all low-lying 
quantum number channels can be determined. The authors use information from 
the overlaps to certain interpolating helds to reliably identify the spin [T3], which is 
challenging due to the rotational symmetry breaking on the lattice. In addition to 
the single-hadron energy levels, the lowest scattering thresholds in the system are 
illustrated by the dashed lines. Above these thresholds, the relation of the observed 
energy levels to hadronic resonances is not straightforward. It should also be pointed 
out that energy levels that arise purely from scattering thresholds can not be distin¬ 
guished from meson excitations when using only single-hadron interpolators. However 
no clear candidates for such energy levels are seen in their simulation, suggesting that 
the number of energy levels observed can be related to mesons with the given quantum 
numbers. 
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Figure 2: Low-lying D (top) and Dg (bottom) me¬ 
son spectra from 


Figure shows low-lying 
positive parity charmed-states 
calculated as part of a study of 
charmed baryons by the QCDSF 
collaboration na. They calcu¬ 
lated the low-lying spectra on 
ensembles with 2-1-1 flavors of 
dynamical quarks (for 2 vol¬ 
umes and 3 pion masses with 
259MeV < mir < 460MeV) and 
the hgure shows results at a lat¬ 
tice spacing of 0.075fm. For the 
D mesons they observe ground 
states corresponding to Du and 
D*7r for = 0+ and 1+ 

which are not shown in the hg¬ 
ure. Their basis was not large 
enough to observe a second Di 

state. For charmed-strange states the authors comment that the 11^1(2536) state 
might be the D*K scattering state. For the with 
ysis uses a only single interpolator and the presence of 
DK scattering threshold can not be excluded. In combination the results pre¬ 
sented here nicely illustrate the pitfalls of using only single-hadron interpola¬ 
tors: The observed overlaps to scattering states seem to depend strongly on 
the technicalities used in the study. This is clearly not a satisfactory situation. 

While in general the rela¬ 
tion of excited state energy lev- m d d. 

els to resonances is not straight- 
forward, Figure [^highlights an- ^ 

other feature from the HSC cal- S - - 

culation of charmed mesons ra- i- L —, ~. « 

Their data clearly reveals hy- . hD- - 

brid meson candidates (contain¬ 
ing both valence quarks and ex- ™ ■ ' - - ~ .■°"- 

cited glue) with negative parity. 

Such states have been conjec- 

,11 • iu 1 Figure 3: Spectrum of D and Dg meson energy lev- 

tured and are seen m the lattice ^ ^ 

1 1 ,. r • . els of negative parity from [121 with hybrid meson 

calculation tor various quantum ^ t- j i j j 

1 energy levels shown in red. 

numbers. 


= O’*" their anal- 
a signal due to the 
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3 Charmed mesons and scattering 


To determine scattering phase shifts and 
the masses and widths of resonances or 
the pole positions of bound states from 
Lattice QCD, Liischer’s finite volume 
method [SI El E! and its various exten¬ 
sions (for a review see [IS]) can be used. 
Within this framework 2-hadron scat¬ 
tering and transitions are well under¬ 
stood and there is some progress for 3 
(or more) hadrons. Just as described in 
the last section, the determination of en¬ 
ergy levels is the first step. Liischer’s 
formula and its extensions then relate 
the energy shifts in finite volume with 
respect to scattering thresholds to the 
phase shift (s) of the continuum scatter¬ 
ing amplitude. 

For charmed mesons this method has 
hrst been applied by taking a look at 
the lowest positive parity D meson reso¬ 
nances in s-wave Dn and D*7r scattering 
PI. The exploratory study on lattices 
with m-jr = 266MeV determined the cou¬ 
plings g rather than T = Figure 

shows the resulting resonance masses 
along with naive energy levels for some 
other D-meson states. Table [T]shows the 
resulting couplings, which are in qualita¬ 
tive agreement with experiment. 

As demonstrated by preliminary re¬ 
sults from the Hadron Spectrum Col¬ 
laboration na, these results can be im¬ 
proved by considering several volumes 
and multiple moving frames. Prelim¬ 
inary results for Dvr-scattering Isospin 
/ = I are shown in Figs. ^and[^ Fig¬ 
ure shows the energies compared to the 
non-interacting energy levels for multiple 
momentum frames. Figure shows the 
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Figure 4: Low lying D meson states from 
[TB] compared to experiment. Only the red 
diamonds are resonance masses, while all 
other states are naive energy levels. 
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D*(2400) 

Di(2430) 

[GeV] 

2.55 ±0.21 

2.01 ±0.15 

gexp [GeV] 

1.92 ±0.14 

2.50 ±0.40 


Table 1: Results for the couplings g from 

PI 
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Figure 5: Preliminary results for the Dvr 
Isospin / = I phase shift from PI consid¬ 
ering the lowest partial wave only. 
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Figure 6: Preliminary results for energy levels of Dvr-scattering Isospin I = \ from 
the Hadron Spectrum Collaboration irzi. 


extracted phase shifts considering only the lowest partial wave for this preliminary 
analysis. A dense coverage of interesting energy region is obtained. 



Figure 7: Left: Energies and pole positions for the D*q( 2317) bound state from [18] at 
two different pion masses compared to experiment. Right: Comparison of the s-wave 
scattering length to predictions from an indirect calculation [20]. 

For charmed-strange p-wave mesons recent results [ISl [TO] including the DK and 
D*K scattering states explicitly into the basis along with quark-antiquark interpo- 
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lators for determinations of the D*g(2317) and i7si(2460) resolve a long standing 
puzzle. A much better quality of the ground state plateau is observed with a com¬ 
bined basis and unlike the quark model predictions and older lattice results the results 
agree qualitatively with experiment (see Figure]^. The resulting S-wave scattering 
length agrees with an indirect determination from EFT-calculations with low-energy 
constants determined by lattice simulations in different scattering channels m , as 
displayed in Figure 

In [21] a reanalysis of the lattice data from [THl HH] is performed by using an 
auxiliary potential V and extracting the parameters of V from the lattice data. 

V = a + (3(s — Sth) T = -^ , 

^ ’ E-i - G 



A generalization of Weinberg’s compositeness condition then yields the probability of 
Ending DK {D*K) components in the physical bound states 

P{KD) = 0.76(12) for the F>:o(2317) , 

P{KD*) = 0.53(17) for the F>,i(2460) . 


The authors further suggest a study with 
increased precision which should also in¬ 
clude the scattering channel. As a 

word of caution, it should also be men¬ 
tioned that the lattice data used suffers 
from sizable discretization effects, so any 
results about the composition derived from 
the data has to be taken with a grain of 
salt. 

Figure shows a prediction of the re¬ 
lated i = \ states with = O^, 1+ in 
the bottom-strange meson spectrum |22j . 
These states have not yet been observed in 
experiment. 

Further interesting results for the scat¬ 
tering of Kaons with D-mesons can be ex¬ 
pected from the Hadron Spectrum Collab¬ 
oration. Figure]^ shows preliminary results 
for DK s-wave scattering in Isospin 1 mi, 
considering only the lowest partial wave. 
Similar results for quantum numbers with 
resonances can be expected in the future. 


Ensemble (2) m = 156 MeV 
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Figure 8: Prediction of the Bs spec¬ 
trum analogues of the Z1 *q( 2317) and the 
77*]^ (2460) from [22] • The uncertainty 
estimate for the bound state poles is 
shown by adding statistical and system¬ 
atic uncertainties in quadrature. 
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Figure 9: Preliminary results for the DK s-wave scattering phase shift in Isospin 
1 = 1 from HU considering the lowest partial wave only. 



Figure 10: Electromagnetic form factors of D and D* as a function of from [23]. 

4 Calculations of qdd^t: and goDp, gD*D*p 

Beyond masses, there is also progress in calculating further properties of charmed 
mesons. A recent example is the calculation of the DD*7t and DDp, D*D*p couplings 
from (transition) matrix elements 

< D{p)\Af,{q)\D*{p,s) > with = u'y^'jid , 

2 2 1 - 

< D{p)\V^{q)\D{p) > with = -cy^c + - -dy^d . 

One recent calculation [23] uses 2+1 flavor gauge conhgurations with a = 0.0907(13) 
and G (300,410,570,700) MeV to extract goD^w, goDp, gD*D*p, the electromag- 
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netic form factors and charge radii of D and D* mesons. Figure W shows the elec¬ 
tromagnetic form factor as a function of the four-momentum transfer Q^. Using a 
vector meson dominance approach with 


MQ") 


9D(*)D(*)n 

ml + Q2 gp 


the couplings 


Qddp — 4.84(34) , 
gD*D*p = 5.94(56) 


are obtained. Notice that disconnected dia¬ 
grams are neglected in this determination. 

Another recent calculation [21] uses 2 fla¬ 
vor gauge conhgurations at 4 lattice spac- 
ings with m,r G (280, 500)MeV to determine 
gDD*TT at the physical point. The authors 
use nonperturbative renormalization and as¬ 
sess the systematic uncertainties in the chiral 
and continuum extrapolations. Rather than 
quoting the result for goD-n directly, the au¬ 
thors report gc which is given by 


gDD*n — -7- gc 

Jn 


Figure 11 shows a compilation of recent lat¬ 
tice results along with the experimental val¬ 
ues and results from QCD sum rules. The 
value from [21] (blue symbol and bar in the 
hgure) leads to T{D*^ —>■ = 50 ± 5 ± 

6keV. 
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Figure 11: Value of gc, comparing ex¬ 
periment, Lattice QCD and Sum Rule 
calculations (from [21] )• The blue 
square is the result from [21] and the 
red circle from [23] . 


5 Searches for exotic charmed states 

A recent study by the HALQCD collaboration investigates doubly-charmed and 
charmed-strange tetraquarks [2S] on 2-1-1 flavor gauge conhgurations with a = 0.0907(13) 
fm and m-,, = 410,570,700 MeV, using the HALQCD method [27]. The search fo¬ 
cuses on bound states or resonances in DD, KD, DD* and KD* interactions with 
havor structure ccud and csud. These contain no quark line diagrams with quark 
annihilation. For the charm quarks a relativistic heavy quark action (a variant of the 
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Figure 12: Left: s-wave 1 = 0 scattering phase shifts in the DD* channel. Right: 
Pion mass dependence of the scattering lengths. Both plots from [25] . 


Fermilab method) is used. The HALQCD method consists of calculating a potential 
as a function of distance r and solving the Schrodinger equation with given V (r) to 
determine scattering phase shifts. An example for the resulting phase shifts is shown 
in the left pane of Figure [T^ The authors obtain a repulsive interaction in all / = 1 
channels and an attractive interaction in all J = 0 channels considered. At the sim¬ 
ulated mir, no bound states or resonances are observed, although the authors state 
that the attraction becomes more prominent at light pion masses and that there is 
some indication that BB* with I{J^) = 0(1’'“) is bound 
In [SB] a preliminary search 


for doubly charmed tetraquarks 
is presented. Results are based 
on a 2 flavor simulation with 
a = 0.075fm and = 490MeV 
and lighter than physical charm 
quark mass. Tetraquarks with fla¬ 
vor structure [cc] [ud] and quantum 
numbers I{J^) = 0(1+),1(1+) 

are considered using a basis of 
tetraquark and meson-meson in¬ 
terpolators, including smeared in¬ 
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Lattice Data i—8—i 

DD* 1.166(4) - 

D(1)D*(-1) 1.230(4) 

-D‘D* 1.254(7) - 

D*(1)D*(-1) 1.314(6) - 

D(2)D*(-2) 1.407(3) 

D*(2)D*(-2) 1.480(6) - 

D(3)D*(-3) 1.659(3) 


terpolators. Figure IR shows some 
of the results. Beyond meson- 
meson states, no additional low- 
lying energy level is observed in 
this preliminary study. 

Further studies search for tetraquarks in systems of two heavy and two light quarks 
with static heavy quarks [281 [291130]. Figure 


Figure 13: Energy spectrum for quantum num¬ 
bers I{J^) = 0(1’'“) from the preliminary results 


m 


14 


shows results from a search for udbb 
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qq = (ud-du)/v'2 qq= uu, (ud+du)/v'2, dd 


qq = (c(1)c(2)-c(1)c(2))/V2 




d in fm d in fm 


Figure 14: Left: Binding energy isolines for two static and two light (u/d) quarks. 
Right: Same for two static and two charm quarks. Plots from [5U] . 


ssbb and ccbb tetraquark bound states |30] obtained by studying potentials of two 
static antiquarks in the presence of two hnite mass quarks. Lattices with a = 0.079, 
a = 0.042fm and ~ 350MeV were used. The ht function used for the Lattice 
QCD potentials is given by 

r(r) = -%xp(-Q^)+n. 


In the system of two static and two light quarks a bound state is observed for the 


scalar channel, as shown in the left pane of Figure 14 


6 Summary and conclusions 

Masses of ground states are well determined in current lattice simulations and 
recent results focus on their properties. For studies of excited states, a large number 
of energy levels can be extracted, including energy levels with a large overlap to 
hybrid meson interpolators [12], providing further evidence for the existence of gluonic 
excitations. Until now very few simulations study close to threshold bound states and 
resonances, but hrst studies [ISl HH] suggest that the closeness of the DK and D*K 
threshold is indeed very important for the Z1*q(2317) and Dsi{2AQQ) mesons, which 
was suggested a long time ago EH as a mechanism for explaining their unexpected 
mass and narrow width. Regarding scattering of D-mesons, preliminary results from 
the Hadron Spectrum Collaboration HT] are very promising. Furthermore, there 
are results from several simulations of explicitly exotic tetraquark states [25l |26l 
I2H1 [2ni EO]- While these searches do not (yet) obtain any evidence for charmed 
tetraquarks, some of the corresponding beauty states might exist [SHlEHlEn]. Beyond 
the current exploratory calculations, an urgent task for the lattice community will 
be better control of (heavy-quark) discretization effects by simulations on multiple 
lattice spacings. 
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